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In this paper we will show that, because of the long-range microscopic memory of the random 
force, acting in the solar core, mainly on the electrons and the protons than on the light and heavy 
ions (or, equally, because of anomalous diffusion of solar core constituents of light mass and of normal 
diffusion of heavy ions), the equilibrium statistical distribution that these particles must obey, is that 
of generalized Boltzmann-Gibbs statistics (or the Tsallis non-extensive statistics), the distribution 
differing very slightly from the usual Maxwellian distribution. Due to the high-energy depleted tail 
of the distribution, the nuclear rates are reduced and, using earlier results on the standard solar 
model neutrino fluxes, calculated by Clayton and collaborators, we can evaluate fluxes in good 
agreement with the experimental data. While proton distribution is only very slightly different from 
Maxwellian there is a little more difference with electron distribution. We can define one central 
electron temperature as a few percent higher than the ion central temperature nearly equal to the 
standard solar model temperature. The difference is related to the different reductions with respect 
to the standard solar model values needed for B and CNO neutrinos and for Be neutrinos. 

PACS number(s): 05.20.-y, 05.30.-d, 73.40.Hm, 71.30.+h, 96.60.K 


I. INTRODUCTION 

After many years of experimental and theoretical work in different laboratories, beyond the demonstration that 
the sun produces energy via nuclear fusion reactions (solar neutrinos have been detected with fluxes and energies 
approximately in the ranges predicted by the known standard solar models (SSM), see Ref.s |l] |[ for the up to date 
state of the art in this field), the accepted conclusions for standard solar neutrinos are the following. 

The solar luminosity is well known and the experimental results (Gallex, Sage, Chlorine, Kamiokande) seem to be 
mutually inconsistent; even if we discharge one of the four experiments. 

Even neglecting these inconsistencies, the measured experimental fluxes of Be and CNO neutrinos are significantly 
smaller than the SSM predictions (we need a reduction factor of about 8 for about 2.4 for <&cno and of 2.28 for 
see Table 1). 

The different reductions of the ‘Be and 8 B neutrino fluxes with respect to the SSM predictions are essentially in 
contradiction with the fact that both 7 Be and 8 B originate from the same parent 7 Be nuclei. Nevertheless, we must 
remember, in relation to what will be stated in the following paragraphs, that in the reaction 7 Be + e~ —> 7 Li + v e 
beryllium nuclei react with the electrons; while in the reaction 7 Be + p —> 8 B + 7, 8 B —> 2a + e + + v e , ‘Be nuclei 
react with the protons. 

The above results could suggest that the hypothesis that nothing happens to the neutrinos after they are created 
in the interior of sun is incorrect. The alternative is that some of the experiments are wrong, this must be checked. 

There are accepted strategies to construct non-standard solar models in order to calculate reduced 7 Be and 8 B 
neutrino fluxes. They are as follows: 

a) To produce models with a central temperature T c smaller than T C SSM . 

b) To play with the nuclear cross sections determining the branches of the fusion chain. However, we must say that 
the resonance 8 He— 8 He has not been found in the very recent measurements of the cross section by the LUNA group 
at Gran Sasso Laboratory Ji| and, as it has been recently shown by Oberhiimmer [?]], we are incapable of solving, in 
the pure nuclear physics framework, the solar neutrino problem. 

c) To take into account correctly the effects of screening the nuclear charges by the stellar plasma (Q. In order to 
understand its contribution to the fluxes and make a comparison with the results from the laboratory experiments 
a new theory of energy loss at very low energies is needed. Very recently Brown and Sawyer have calculated the 
screening effects in stellar plasma and have found that the enhancement factor to the fusion cross sections is very 
close to one [()]. 

d) To calculate from the beginning the opacity, taking into account plasma collective effects, as done recently by 
Tsitovitch and collaborators @ and by Ricci |TT|]. However, we must remember that, in this case, the reaction rates 
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must be recalculated if the electron distribution deviates from the Maxwellian one due to the presence of particular 
collective effects, external forces or other factors. 

Generally it is accepted that for standard neutrinos the actual experimental results seem to be in disagreement 
with SSM and with all the solar models we could build even considering neutrino oscillations Ell¬ 
in spite of these arguments, we wish to discuss, within the standard solar model description, the thermostatistics 
to which ions and electrons of the solar core should adhere. In fact, we want to show that the Maxwell-Boltzmann 
statistics must be substituted with the non-extensive statistics, suitable for particles subjected to long-range micro¬ 
scopic memory. 

The resulting statistical distribution is, for the solar core, a very modest modification of the standard Maxwell- 
Boltzmann distribution p~i| . |i~5| |. 

The memory acts on the dynamics of the particles of the solar core because the electron mass is much smaller than 
the mass of the other constituents of the core plasma and because the proton mass is smaller than that of light and 
heavier ions. The normal memory is related to a time correlation of the random force which is a delta function. In this 
case the diffusion is normal, the equilibrium distribution is Maxwellian. The long-range memory implies a random 
force with time correlation different from a delta function as suggested by the results of measurements of neutrino 
fluxes. 

The diffusion coefficient of a system of particles, subjected to these forces, is anomalous (non-Brownian) and the 
equilibrium distribution function is a non-extensive Tsallis distribution [ P~6 |JTt[ . 

We describe, in the following sections, the main features of the non-extensive Tsallis statistics, the normal and 
anomalous diffusion coefficient, a few applications of astrophysical interest, the application to the solar neutrino 
production using the results from the SSM of Clayton and collaborators |l8|]ig| l and give the physical motivations of 
the validity using, in this case, the non-extensive statistics. 

In the Appendix A we discuss the anomalous diffusion coefficient in the framework of the Tsallis statistics because 
this quantity represents the link between the distribution adopted and the motivation of its use (i.e. the long-range 
time correlation of the random force). 


II. THE USE OF NON-MAXWELLIAN DISTRIBUTIONS IN THE RECENT PAST 

A customary mode of procedure to change the SSM nuclear rates can be based on the use of the appropriate 
statistics (the generalization of the Boltzmann-Gibbs statistics or the Tsallis non-extensive thermostatistics 
which describes the behavior of the different particles composing the solar core. In this work we introduce a modified 
Maxwellian function which describes an equilibrium distribution and is the first order approximation of the Tsallis 
distribution function. 

The physical motivations of this choice are of a many-body nature as it will be explained later. This choice allows a 
reduction of the Maxwellian rates and a progression towards calculations of the neutrino fluxes not in conflict with 
experimental results. 

To change the statistics and, of course, the distribution function means also to change the energy spectra of 
the neutrinos produced. Possibly, in future experiments at Gran Sasso Laboratory |f2(j| or, for instance, in Hellaz 
experiment fjlj i spectra could be measured with enough precision to allow the verification of the validity of the 
generalized, non-extensive statistics in the solar plasma. 

In the past, after the suggestion of Kocharov |^2|, Clayton and collaborators E-El calculated the evolutionary 
sequences for the sun under the assumption that the relative energies of the ions are not exactly Maxwellian. 
Depletion of the high energy tail of the distribution function of relative energies in excess of 10 kT produces a marked 
reduction in the counting rate of the Cl experiment, while producing minimal changes in the usual solar models. 
Those changes that do occur can be compensated for by small changes in the initial helium concentration. 

The mechanism leaves its own signature in the nuclear rates and in the energy spectra of solar neutrinos. 

Clayton could not calculate a physical cause of the deficiency in the number of energetic pairs and suspected that 
the problem lies in many-body physics: long range Coulomb interaction conspires in some unknown way to quench 
relative energies above « 15 kT. 

Correction to the Maxwellian distribution function is given by the factor exp[— 6(E/kT) 2 ]. Clayton found for the 
solar neutrino problem S = 0.01 (E is the c.m. energy). 

Departure of the two-particle distribution from the Maxwellian distribution function in the solar core is not un¬ 
physical. It should be possible and, as we will see in this work, the reason is long-range microscopic memory E,^4j , 
more than the long-range interactions (Coulomb and gravitational). 
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Haubold and Mathai [p5| , ^6| , avoiding mathematical approximations, derived closed form representations of nuclear 
reaction rates including resonant, non resonant, screened cases in terms of special functions. All these calculations 
have been accomplished without attributing an exact and precise physical meaning of the depletion of the Maxwellian 
distribution function used, except the following, that many-body effects in the solar core should be taken into account. 
We also performed in the recent past calculations of nuclear rates and recombination cross sections for opacity using 
non-Maxwellian distributions j2?],|28|], before our first application of non-extensive statistics to solar neutrino probem 


The inclusion of the helium diffusion in solar models has been widely studied by Bahcall and collaborators in the 
frame of Maxwell-Boltzmann statistics J29|,|30|]. 

In this paper, because of physical motivations, we want to introduce, for the description of the solar plasma, a 
generalization of the Boltzmann-Gibbs statistics and by means of a reevaluation of the nuclear rates of interest to 
calculate the neutrino fluxes. 


III. NON-EXTENSIVE STATISTICS 


Recently, non-extensive statistics (generalization of Boltzmann-Gibbs statistics), based on a new generalized defi¬ 
nition of entropy, becoming the well known standard definition when the characteristic parameter q —> 1 , is widely 
considered in many different physical phenomena. 

The distribution function of these statistics has not only a depleted high energy tail, but in reality vanishes when: 

kT 

E = — if q < 1 • (1) 

Tsallis Jl(| introduced in 1988 this generalization of thermodynamics and statistical physics to describe systems with 
long-range interaction (gravitational) and long-range memory [^l] [l3j . 


Let us describe briefly Tsallis statistics. 

For a system with W microscopic state probabilities pi > 0, normalized as Pi = lj we have 2 axioms: 

1 ) the entropy is 

w k 

S q = J2^ 5,* = —P?" ) > (2) 

»=i ^ 

for q -> 1, Si = -kYjiPdogpi. 

2 ) given an observable O with Oi eigenvalues, the mean value is 

O q = EZiP>i, q^l: O^TZiPiOi- 

The validity of the two axioms lies in the comparison with experiments. 

Properties are the following: systems with q < 1 give more weight to rare events, systems with q > 1 give more weight 
to frequent events. 

The generalized entropy S q is positive, the microcanonical ensemble has equiprobability, entropy is concave (q > 0), 
convex (q < 0); Legendre transformation structure of thermodynamics is invariant for all q. 

Additive rule (non-extensivity) is : 

S q (AuB) = S q (A) + S q (B) + {l-q)S q (A)S q (B) . (3) 


The Tsallis distribution function is given by the expression 



1 - (1 - «) 


Ei_ 

kT 




and 
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( 4 ) 

( 5 ) 


Application of generalized non-extensive statistics are the following subjects: matter distribution of self gravitating 
systems, turbulence, anomalous diffusion, cosmological models, big bang background radiation We describe 

here only those applications related to astrophysical problems (see Sections 5, 6 and 7). 

In Fig.l the behavior of the function pi is shown for values of q greater and smaller than one. 
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IV. NORMAL AND ANOMALOUS DIFFUSION COEFFICIENTS 


Is is well known that neutrinos of the four solar fluxes <l>p, $B e , &cno and are produced in an equilibrium 
plasma in the solar core. In fact, the average time T nuc i between two fusion reactions is much greater than the average 
time between two Coulomb collisions: T nuc i Tcouiomb- We can define 


1 1 

7~nucl — “ “ i 

Tin < O V > nucl 


( 6 ) 


where n n is the nuclear density and < a v > nuc i is the nuclear fusion rate averaged over a Maxwellian distribution. 
The quantity Tcouiomb is related to the Brownian diffusion coefficient to] 


n _ kT 

TJ Br — ^Coulomb 

m 


( 7 ) 


A factor of 10 20 is in favor of equilibrium and Maxwell-Boltzmann distribution function is used to describe both 
electrons and ions dynamics and reaction rates. In support of the validity of the Maxwell-Boltzmann distribution, 
explicit calculations in which the Fokker-Planck and Boltzmann equations were solved show that the Maxwellian tail 
is filled, after few Coulomb collisions, when the different colliding particles (both electrons and ions) are all considered 
Brownian. 

It is known that neutrino production nuclear rates, reduced in respect of their SSM values, might lead to an 
evaluation of the solar neutrino fluxes in agreement with the solar luminosity and not too far from the results of 
measurements of fluxes accomplished in underground experiments. 

Distribution functions of equilibrium are not only Maxwellian distributions. If the diffusion coefficient is not 
Brownian (or normal) because the medium is dense enough (like the solar core) and the random force has a long-time 
correlation, then the equilibrium steady state is described by a distribution which is not Maxwellian, although it may 
be very close to a Maxwellian. 

Recently, we have shown that, in the energy space, we must introduce an anomalous diffusion coefficient to solve a 
generalized Fokker-Planck equation, given by the expression 


D = D B r 
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( 8 ) 


valid up to the energy E = kT/{l — q)\ to this coefficient we can relate the anomalous collision time 
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P <t (TV > Coul 




E ' 
kT 


( 9 ) 


where p is the average density. 

We are still in equilibrium conditions, but the distribution is a depleted Maxwellian and of non-extensive Tsallis type. 

If the ions (proton, helium, light and heavy-ions) and the electrons can be considered Brownian particles, their 
diffusion coefficients are normal and their stationary distributions are Maxwellian. This description is correct also in 
the presence, in the dynamical equations, of a random force with a time correlation equal to a delta function [^| ^6| . 
The use of the distribution appropriate to the non-extensive (q ^ 1) thermostatistics has the physical meaning that 
the description of the dynamics of the solar core constituents contains anomalous, rather than normal, diffusion co¬ 
efficients. 


In the solar core we must add, in the dynamical equations of the different constituents, a random force F(t) with 
time correlation different from the delta function (see the Appendix A for the needed details): 


<F(0)F(t) >=F 0 {(3)t~P 


( 10 ) 


where /3 is a parameter, in this case, smaller than one and related to q by 


2 -q 1 + 2 S 


( 11 ) 


(for (3 = 1, or q = 1 and 5 = 0 the delta function time correlation is recovered). 

The choice of a non-delta time correlation is suggested by the comparison of calculations with the experimental results, 


4 










i.e. by the values of the parameter q required by the results of the measurements and, consequently, by the values of 
the parameter (3. 

In the solar core, the anomaly consist in a subdiffusion, because the value of the parameter q to be used lies in the 
range: q = 0.952 4- 0.994 (or S = 0.003 4- 0.024), as explained in the following Section 7. 

The subdiffusion correction to the normal behavior is very small but sufficient to calculate neutrino fluxes in agreement 
with the measured fluxes. 

In the solar core, the anomaly is not due to the long range gravitational interaction, rather it is due to the long-range 
memory of the random force. 

In the Appendix A we discuss the meaning of a non-linear or non-Brownian diffusion coefficient and the solutions 
of the Fokker-Planck equation appropriate to the non-extensive thermostatistics. 


V. STELLAR POLYTROPES 


This is the first example of application of non-extensive statistics. 
The equation of state appropriate to the stellar polytropes is 

V = KLp< , 


( 12 ) 


where V is the pressure, p the density, 7 a constant related to specific heats defined by 7 = 1 + 1/n (n is the polytropic 
index). 

The idrostatic equilibrium equation has spherical symmetric solutions corresponding to compact spherical configura¬ 
tions of self-gravitating mass (stellar polytropes). 

Relation for stellar polytropes between n and the Tsallis parameter q is |M 


3 

H ~ 2 


1-9 


(13) 


The index n must exceed 1/2 to avoid singularity in the gamma function; n > 5 gives rise to unnormalizable mass 
distribution (unphysical). The Tsallis distribution function has a spatial cutoff of the mass distribution (compact 
nature of the stellar polytrope). 

In the solar core we have a constraint which fixes differently the relation between n and q |l-l]Jl5t| 


9 = 1- 


n + 1 


T = 


kTp 

~v~'- 


q=l-25 


(14) 


the quantity r is very close to one if an ideal gas of particles having a certain average molecular weight is considered; 
in fact, the equation of state of an ideal gas, within the Tsallis statistics, differs from the classical one for the addition 
of terms whose contribution depends on powers of S (in the limit q —» 1 , this contribution can safely be disregarded) 


VI. APPLICATION TO GALAXY CLUSTER VELOCITY 


We have shown that the observational data recently provided by Giovannelli et al. J48| (COBE) and discussed by 
Neta Bahcall and Peng Oh 1491 and by Moscardini et al. |30| concerning the velocity distribution of clusters of galaxies 
can be naturally fitted by a statistical distribution which generalizes the Maxwell-Boltzmann one IsTJ]. 


In this generalization of Boltzmann Gibbs statistics the probability function of having a cluster with velocity greater 
than v can be written as 


P(> v) = 


[ dv [1 - (1 - q)(v/v 0 ) 2 ] 

J V 


21 ?/(! — 9 ) 


dv [l - (1 - q)(v/v 0 ) 


21 <?/(!—<?) 


'Umax — 


^0(1 - 9) 


-a I ” 1 / 2 


if q < 1 
if q > 1 . 


A remarkably good fitting with the data is obtained for q = 0.23±q q 5 and vq = 490 ± 5 km/s. 
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VII. APPLICATION TO SOLAR NEUTRINO PROBLEM 


Taking advantage of the calculations accomplished by Clayton and collaborators, based on their solar model (stan¬ 
dard) which uses a distribution function of Maxwellian type with a depleted tail (this is a first and suitable approx¬ 
imation of the Tsallis distribution), we discuss the application of the non-extensive statistics to the solar neutrino 
problem. In Refs @0 we report how the constraints imposed by the solar core compel us to derive the relation 
(14) instead of Eq.(13) suitable for stellar polytropes; results of the application of the Tsallis statistics to solar neu¬ 
trino problem are shown without mentioning the physical motivations to use this new statistics (mainly due to the 
anomalous diffusion of light particles rather than normal diffusion appropriate to heavier ions). 

A value of the Clayton parameter 6 different from zero (we remember that <7 = 1 — 26) makes the star more luminous 
and reduces the rate of energy production at a given temperature. The solar core contracts to higher temperature. 
However, the increase of temperature at given solar luminosity does not increase the neutrino fluxes that decrease 
with 6 (except <h p , but its growth is very slow and its value is almost constant) because of the depleted Maxwellian 
tail and the consequent rearrangement of the particle distribution. 

We introduce the reduction factor f red of the B, CNO and Be fluxes as a function of 6 (or q) 


^B,CNO,Be 



SSM 

B,CNO,Be > 


(15) 


and the enhancement factor fi ncr of the proton flux as a function of 6 (or q ) 


% = 


fir, 


(fjSSM 

P 


(16) 


In Fig. 2 we show both the reduction factor f red and the increasing factor fi ncr as function of 6. 

These two quantities are derived from the calculation of the neutrino fluxes as function of 6 by Clayton and collabo¬ 
rators. The SSM fluxes are all normalized to one at 6 = 0. 

We have obtained an expression for the factors valid up to <5 = 0.025 (an extrapolation of the results reported by 
Clayton and collaborators, in the range between <5 = 0.020 and <5 = 0.025, has been done): 


fred{B) = 1 + 92581.75(5 1 92 , 
f red (CNO) = 1 + 23.6(5° 49 + 305 , 

fred(Be) = 1 + 851.7L5 1 - 29 , (17) 

fincrip) = 1 T 0.89(5°' 92 . 


We can clearly see that the selection of 6 = 0.003 (q = 0.994) produces the correct reduction and enhancement 
factors (see Section 1) of the SSM values of &cno and <f> p and then neutrino fluxes in agreement with the 
measured values and with the solar luminosity; the choice 6 = 0.024 (q = 0.952) is the right selection for 
The reduction factors f re d have been derived in such a way that the four fluxes satisfy the constraints imposed by 
the solar luminosity and by the results of measurements of the Gallex, Sage, Chlorine and Kamiokande experiments, 
within their precision |^,[|]. 

The need for two different choices of 6 is due to the Be flux and depends on the electron distribution, while the B 
flux (and, of course, the proton and CNO fluxes) depends on the proton distribution. 

In fact: 

oc N{e~)N(f Be) < e~ 7 Be > and oc N(p)N( 7 Be) < p 7 Be >, 
where N(e~) is the electron density, N(p) the proton density, < e~ 7 Be > the nuclear capture rate electron-beryllium 
and < p 7 Be > the nuclear rate proton- beryllium (analogously we can say of the other two fluxes). See the Table 
where the results on the fluxes are reported. 

The consequences of this derivation is that, in the solar core, electrons behave differently from protons and from the 
other light and heavy ions. The latter are nearly Maxwellian, the electron distribution must have a more depleted 
tail. In addition, the polytropic index n of electrons differs from that of the other constituents. The solar core 
plasma model, consisting of two main components (electrons and ions), with different features, is under study. Here 
we still want to examine the effect of these physical facts on the central solar temperature by using the well known 
dependencies of the fluxes on central temperature as given in the literature [|^p9p(i|]. 
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VIII. TEMPERATURE OF THE CENTRAL CORE 


We derive the dependence of the four fluxes upon central temperature T c when all SSM calculations are executed at 
T C SSM , taking advantage of the expressions reported in Ref. |4||3C|. We have derived from Clayton and collaborators 
calculations the expression 


T c = T < f SM (l + 3.125) , 


(18) 


valid for S < 0.1. 

The dependence is as follows 

oc (1 - 0.88f n )(l + 0.44f°' 62 ) , 
oc T 25 [l + 10 4 (t — l) 1 ' 92 ] -1 , 

$cnooc(xT™+ (l^x)T™)[l + 13M{t-l) 0A3 }- 1 , (19) 

^BeOcT^fl + lQe^-l) 1 - 29 ]- 1 , 

where t = T C /T C 5SM = 1 + 3.12 5, x represents the percentage of the contribution to the CNO flux from the ions 13 N 
and 15 O and 1 — x the contribution from 17 F. 

We can see that T c increases by 1 per cent in respect to T^ SM using the value of <5 appropriate to ® p ,b,cno fluxes 
(this temperature is the ionic plasma temperature tJ^ = 1.01 T C SSM ), and increases by seven per cent using the value 
of 5 appropriate to &Be (this temperature is the electron plasma temperature Tc^ = 1.07 T^ SM ). 

The meaning of the results of the above relations is different from that of the results coming from the relations of 
other authors. While, according to these authors, the central temperature must decrease in respect to the SSM value 
to reproduce the behavior of the ratio of the fluxes, in our approach the central temperature T c must slightly increase 
with S and the fluxes (except which is almost constant) decrease because of the depleted tail and the consequent 
rearrangement of the particles. 

Helioseismology is probing the interior structure and dynamics of the sun with great precision |5^-|54|. Acoustic 
waves and internal gravity waves are strongly influenced by the central temperature and by the structure of central 
regions. The reliability of the non-extensive statistics applied to the solar core can be tested in the context of 
helioseismological studies. This is the topic of our investigation in the near future. 


IX. CONCLUSION 

The validity of the description of the solar core shown in this work, by means of a generalized Boltzmann-Gibbs 
statistics, could be verified in future planned experiments by measuring the electron central temperature T c , through 
the detection of the Be line and by measuring the neutrino energy spectra, looking for a shift of the Gamow peak and 
of the neutrino maximum energy due to non-extensive thermostatistics effects, if very high precision will be obtained. 
It could be of practical use, from a theoretical point of view, to evaluate the effect of the neutrino oscillations of the 
fluxes in the context of the non-extensive statistical description of the solar core. 

We have shown that anomalous diffusion (different from Brownian) of electrons and of protons, caused by the long- 
range microscopic memory of the random force in the solar core plasma, is the cause of the stationary distributions 
of the non-extensive statistical type we must use. Light and heavy ions have a normal diffusive behavior. As a 
consequence, the values of the nuclear rates of the different reactions decrease in respect to Maxwellian evaluations 
(we have taken advantage of the SSM of Clayton and collaborators). 

A particular choice of the value of the non-extensive Tsallis parameter q for the electrons (q = 0.952) and for the 
protons (q = 0.994) allows the calculation of neutrino fluxes. The values of q differ from the value q = 1 (Maxwellian 
distribution) for very small magnitudes. The value q = 1 is assumed for light and heavy ions. 

This description implies that the electrons have a behavior which differs from the Maxwellian behavior more than that 
of the protons, while we can assume that light and heavy ions have a normal behavior, or normal diffusion coefficient, 
which implies Maxwellian distribution. 

Our results, based on the standard solar model developed by Clayton and collaborators, imply a weak increase of the 
central temperatures (both electronic and ionic), the fluxes decrease because of the depleted tail of the distribution. 

Finally, we must remember that the validity of the non-extensive statistics in the solar core plasma is also based on 
the polytropic nature of the sun, whose polytropic index n has a finite value with the consequence that the distribution 
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must generally be non-Maxwellian. 

The approach described in this work will be refined in order to take into account the information given by the 
helioseismological studies and measurements. 

Our approach is based on a standard solar model with diffusion (anomalous and normal) of the constituents of the 
solar core plasma, therefore it is highly possible that the test with the helioseismological data be positive, but it must 
be checked; this will be accomplished soon. 

We realize that the results reported in this paper are not definitive. In fact we do not indicate errors and precisions 
of the figures we have derived and shown. This is due to the reason that we do not have used a solar model code 
and we have not accomplished calculations abinitio. However, we have taken advantage of a standard solar model 
previsions and our results are an indication of a possible solution of the solar neutrino problem (anomalous diffusion 
for the lighter constituents of the core). We hope that in the next future the non-extensive Tsallis statistics will be 
tested within one of the standard solar models actually in use. 


X. APPENDIX A 

In the solar core we can distinguish the following constituents: 
protons; helium ions; other light ions and heavy ions and electrons (with TOj on mue > Wp rn e ). 

Due to the magnitude of m e , the electrons diffusive behavior is a little less normal than that of light ions, because their 
dynamics, described for instance by a generalized Langevin equation, require the inclusion of a particular random 
force. For the same reason, the magnitude of their mass, heavy ions have a normal diffusive behavior. We assume 
that also helium has a normal diffusion. 

More explicitly, the deterministic equation mv + av = 0, whose solution is v(t ) = c(0) exp(—yf) (7 = a/m = t _ 1 ), is 
not sufficient to describe the diffusion of particles with a mass smaller than the masses of the particles composing the 
rest of the plasma (their velocity due to thermal fluctuations is considerable). 

If the time correlation of the random force F(t), which we must add, is a delta function, the equilibrium distribution 
is also a Maxwellian distribution. 

From another aspect, a long-range microscopic memory, i.e. a time correlation of the random force different from a 
delta function, produces an equilibrium distribution of a non-extensive statistical category. 

Subdiffusion is responsible for the depletion of the Maxwellian tail. The equilibrium distribution of electrons and 
of protons and helium ions can be derived from dynamical equations containing a colored noise, more coloured for 
electrons than for proton and for helium ions. 

Anomalous diffusion, like subdiffusion, is related to the normalized equilibrium velocity autocorrelation function 
C v (t) @ 

Cv(t) = * < v(0)v(t) > , ( 20 ) 

<v z > 

where v(t) is the velocity of the particle at time t and C v {t ) does non depend on the time origin. 

Mean square displacement is 


< x(t ) 2 >=t f C v {r)dT — f tC v (t)cIt 
Jo Jo 


( 21 ) 


If C v (u) decays faster than u~ 2 after a long period of time, one observes normal (or Brownian) diffusion. 

We are interested here in subdiffusion because our analysis of the measured neutrino fluxes suggests that this is the 
type of diffusion of the solar core particles. This regime is reached when C v (u) ~ Au~^ 2 ~^ as u —> 00 (0 < (3 < 1, 
A < 0 and / 0 °° C v (u)du = 0) then < x(t ) 2 >« Bt 13 as t —> 00 . 

If the particle is moving in the positive x direction, it is more likely to move in the negative x direction in the next 
instant. The fluctuating velocity reverses its direction very often. 

The time correlation of the random force can be written, for instance, as 


<F{0)F(t) >= jex P (-t/d) . 

If 7 ? —> 0 the correlation is given by 2D5(t) which refers to a normal diffusion (therefore q is linked to d). 
Following Wang [M , we can write generally 


( 22 ) 
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( 23 ) 


<F(0)F(t) >= C f (t) = F 0 (P)t-f> 


We relate now the parameter f3 to the Tsallis parameter q. Tsallis and Bukmann |3(| have shown that correlated 
anomalous diffusions are described by the following generalized Fokker-Planck equation in the space coordinate p6| 


= ~^{F{x)\jp{x,t)Y} + D-^ \p(x,t)] v 


(24) 


where F(x) = —dV/dx is an external (drift) force associated with the potential V(x) (F(x) = k\ — IC 2 X , /c 2 > 0). 
The Tsallis parameter q is related to p and v 


q = 1 + p - v 


(25) 


As here we want to preserve the norm of the distribution, we must pose p = 1. 

Therefore: q = 2 — i/ = 1 — 28 and then v = 1 + 28 (8 is the Clayton parameter). By comparison, we have the relation 


0 = 



1 _ 1 
1 + 2 S~v 


(26) 


The normal diffusion is recovered when (3 = 1 (y = 1), i.e., as we know, when q = 1 (8 = 0). 
The Fokker Plank equation for the particles of the solar core is (p = 1, v = 1 — 28) 


^p(x,t) = —^{F(x)p(x,t)}+D-^\p(x,t)) 1+25 . (27) 

The stationary solution of the above equation is the Tsallis distribution in the coordinate space (x has the dimension 
of an energy square root) 


p(x) 


f Af(q, kT) 

1 (1-9+ 2 

kT 

1 o 






otherwise 


where A f is the following normalization constant 


Af{q, kT) = 7T 1//2 



r 


(5-3q\ 



(28) 


(29) 


The results obtained in the Appendix can be recovered using the approach described in Ref. Jd2[ where a new 
definition of the diffusion coefficient in the momentum space, whose expression is reported in Eq.(8) of the text, was 
introduced. 

In the momentum space the equilibrium distribution can be derived by means of the Fourier transformation of the 
distribution p{x), following the procedure outlined by Tsallis et al. |36|| (its analytic expression is in terms of the 
modified Bessel function 1C (^5|). In the situations where g —► 1 (as in the case of the solar neutrino problem) we can 
easily verify that the Tsallis distribution in the coordinate space, given by Eq.(28), corresponds to another Tsallis 
distribution in the momentum space. In fact, for q —► 1 we can write (here x is adimensional): 


p{x) (x[\ — — q)x 2 } 1 ^ 1 ^ ~ exp(— x 2 — 1 ^ x 4 ) . (30) 

The Fourier transformation of p(x) is the Clayton distribution 

F\p(x)\ ~ exp-T[i + (i -g)^j- - (1 -g)^] ~exp[-^- - (1 -g)^] , (31) 

which approximates, in the momentum space, the Tsallis distribution 

F\p{x)\ tx [1 - (1 - g)p 2 /4] 1/(1 “ 9) ■ (32) 
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Figure Caption 
Fig.l 

The Tsallis distribution function p(Ei) for different values of q: q = 0, 0.5, 1 (MB), 1.5. The typical depletion of the 
high-energy tail of the Maxwellian (q = 1) is clearly evidentiated for q < 1. 


Fig.2 

The reduction factor f re d (B, CNO, Be) and the increasing factor finer (p) as functions of the parameter 5 (q = 1 — 2 6), 
up to S = 0.025. We call the attention to the fact that SSM values of the B , CNO and Be fluxes must be divided 
by fred (B, CNO, Be) to obtain the fluxes within the non-extensive statistics, while the SSM proton flux must be 
multiplied by fi ncr (p)- All SSM fluxes are normalized to one at S = 0. 


Table Caption 
Table 1 

The SSM fluxes (8 = 0) from Ref. | |30|| , the flux = finer® p SM and the fluxes ®B,CNO,Be = ®B S CNO,Be/ fred 
evaluated within the non-extensive statistics. All fluxes are in 10 9 cm' 2 s _1 but >!>b in 10 6 cm' 2 s -1 (we recall that 
the Kamiokande experiment gives the result 4 >b = 2.9 ± 0.4). 


(5 

®p 

•hs 

®CNO 

®Be 

0 

59.24(1± 0.01) 

6.62(1™) 

1.17(1±0.2) 

5.15(1±0.06) 

0.003 

60.66 (f incr = 1.024) 

2.85 (f red = 2.326) 

0.48 (/red=2.46) 


0.024 




0.65 (f red = 7.93) 
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